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Elliptical Solutions to the Standard Cosmology Model 
with Realistic Values of Matter Density 


Ahmet Mecit OztaH Michael L. SmithU 


We have examined a solution to the FRW model of the Einstein and de Sitter 
Universe, often termed the standard model of cosmology, using wide values 
for the normalized cosmological constant (U A ) and spacetime curvature (U^) 
with proposed values of normalized matter density. These solutions were 
evaluated using a combination of the third type of elliptical equations and 
were found to display critical points for redshift z, between 1 and 3, when P A 
is positive. These critical points occur at values for normalized cosmological 
constant higher than those currently thought important, though we find this 
solution interesting because the U A term may increase in dominance as the 
Universe evolves bringing this discontinuity into importance. We also find 
positive U A tends towards attractive at values of z which are commonly 
observed for distant galaxies. 
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1 INTRODUCTION 

The standard view of Universe expansion, the Friedmann-Robertson-Walker 
(FRW) model, seems to mimic our current situation well given the approx¬ 
imation of an isotropic and homogeneous sprinkling of matter across large 
dimensions and inclusion of vacuum energy. The standard model has several 
variants and most include significant positive values of A the vacuum energy 
constant; this energy of deep space seems to encourage Universe expansion 
and act as anti-gravity. Some current detailed models predict a greatly in¬ 
creasing importance of this energy with continued Universe expansion, with 
suggestions we are entering the age of dominance by vacuum energy (Be- 
har and Carmeli, 2000). Another recent, interesting model details such an 
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energy coupled to absolute time, allowing an initial vacuum energy driven 
inflationary phase immediately after release from singularity, followed by a 
more leisurely relaxation of this energy (Bisabr, 2004 and references therein). 
Perhaps the current vacuum energy is a residual of the first moments of the 
Universe. 

Recent measurements of supernovae type la (SNe la) distances and as¬ 
sociated redshifts z (or receding velocities), have uncovered the existence 
of a significant current value for the cosmological constant (Tonry et al., 
2003, Reiss et al., 2004, 1998). The measurements of SNe la distances and 
receding velocities are subject to a myriad of possible experimental errors 
demanding careful data collection and detailed analyses, for evidence of this 
positive repulsive energy is obtained from quite distant SNe la explosions. 
Still these SNe la studies present data with the smallest experimental er¬ 
rors to date and are the most useful for model testing. These authors also 
take the positive A to mean the Universe is expanding more rapidly now 
than in the past and the local Hubble constant, Hq is gradually increasing 
with time. Without this energy from deep space we should expect Hq to 
slowly decrease over time due to the self attraction of matter and energy. 
Knowledge of Hq is also important because the upper bound of the age of 
the Universe is fixed by the local value of Hq in many models. Surprisingly, 
Hq is not known to better than ±10% (Freedman, et al., 2001), indeed the 
Universe age is not known to better than ±10% if the age of low metal stars 
is a gauge (VandenBerg, et al., 2002, Grundahl, et al., 2000). One reason 
for this is the dependence of Hq upon matter density which itself is highly 
dependent upon the distances between stars and galaxies and dust density. 
Another reason for inaccuracy being the relative time span of our data, 
both SNe la and star elemental composition, are instantaneous respective 
to processes of the Universe. 

We have examined the standard model of the expanding Universe using a 
solution of the elliptical form and found that at realistic matter densities 
this solution exhibits discontinuities at values of the normalized cosmologi¬ 
cal constant slightly larger than current estimates. When we examine this 
solution for the smaller normalized matter densities of the future (and val¬ 
ues suggested a few decades ago), this discontinuity presents problems for 
vacuum energy dominated models. Use of the cosmological constant in the 
standard model also clouds predictions of galactic redshifts at epochs only 
slightly older than those collected in SNe la studies. Though it would be 
useful to be able to calculate distant galactic ages from redshift alone, now at 
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7 or higher (Kneib, et al., 2004), such are unreliable if the positive vacuum 
energy is significant. We do not find discontinuities for negative vacuum 
energies (attractive) at realistic values of normalized matter densities. 

2 THE FRIEDMANN MODEL 

We use the conventions of Carroll et al. (1992) with a Friedmann equation 
modeling our expanding Universe 

~ ,R. o 8irG A k 

~ ~!T Pm+ 3" ~~ R?' 

Here the p m is the matter density and the first term attractive, while the 
second term may either be repulsive or attractive, it is usually meant as 
repulsive. The constant of integration k, may take values of -1, 0 or +1, 
for a Universe with ’’open, flat or closed” geometries, respectively. The 
H represents the Hubbell parameter with R the expansion factor for the 
evolving Universe. 

Einstein first proposed his gravitational equation without a cosmological 
constant and preferred a very slightly closed Universe with matter domi¬ 
nating (Einstein, 1915), that is k > 0. He later introduced the cosmologi¬ 
cal constant into his gravitational equation after Friedmann and Lemaitre 
pointed out the Universe, populated with considerable matter, should either 
be expanding or suffer contraction, but astronomers had not yet firmly dis¬ 
covered other ’’Universes” or galaxies outside the Milky Way. This constant 
allowed the possibility of a static Milky Way (and Universe), which was the 
limit of knowledge early in the 20th century. Hubble, Wirtz, Slipher and 
others later pointed out that most other galaxies were following trajecto¬ 
ries away from the Milky Way with the implication that the Universe has 
no possible stationary reference point, in confirmation of Einstein’s propos¬ 
als. It seems that Einstein later regretted introduction of this cosmological 
constant, nonetheless, this concept has recently regained popularity in cos¬ 
mology to explain certain observations. 


3 THEORY 


It is common to introduce normalized parameters for matter, vacuum energy 
and geometry 
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We use the typical conditions of normalization across these three parameters 
following the convention of Carroll et al. (1992) 

1 = Q m + £l A + flfc (3) 


where Q m represents normalized matter density, normalized vacuum en¬ 
ergy density and Q k normalized (possible) spacetime curvature; the radiation 
density term, being small at present, is included with matter density and 
we will briefly review the equations of interest. If we allow p mo to represent 
the present matter density which is ( 47r / 3 ^ 3 ) i n equation (1) to give us the 
FRW model at the present time we have 


dR 

Rdt 


H: 
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Now we substitute for with a(t ) 


d (t) \2 _ tj 2 Qm 

a(t) ~ ° a(t) 3 
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Q k 
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and then multiplying through by "iX 
state 

1 da(t) 2 _ Q m 
H'l dt a(t.) 


gives us one equation of our current 

+ ft(t)"fl a T blfc. (6) 


This allows us to introduce the dimensionless parameter r for time as r = 
H 0 t 

2 = _1 ^ + 3 + ^ ( 7 ) 

dr a(T) 

da( ' T " > = + a(r) 3 + flfca(r). (8) 

dT V a ( T ) 

Then inverting and integrating both sides from the past t\ to the present r G 
we have _ 



y 7 o(r) 


la 1 a(r ) 3 + ^ka{r) 


da{r) 


(9) 


and substituting 1/(1 + z) for a we have the integral from the past z\ to the 
present of 0 we arrive at 


T 0 - n = - 



1 

(1 + z)y 7 n m (l + z) 3 + Ha + blfc(l + z) 2 


dz. 
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This is similar to an equation modeling redshifts as presented in Peebles 
(1993). We shall use the variable y for 1 + z and bring out * from within 
the integral of (10), so the integral becomes 


t 0 - n = 


y/^m ■ 


' 1+21 


y 3 + 




+ 


n k y 2 


--dy. 


( 11 ) 


We shall now change the variable once again, allowing y = u — -^-H m which 
also changes the integration limits of the following equation 


" 1+21 + 
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and by using a similar substitution we can simplify the denominator in terms 
of u with the same limits of integration as equation (12) 

To 

The right side of equation (13) is solved as an elliptical function of the first 
and third type, where 


Tl = 


y/^rn j (U - 3 rr)y/{u - Ui)(u - u 2 )(u - u 3 ) 
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To ~ Tl = 


y/u:i - ui(u 3 - 


x n(x, 12, k ) 


(15) 


and the function can be evaluated using the limits of Eq. (12). We may also 
proceed from equation (10) by inverting the limits of integration, changing 
sign and substituting for r 0 — r\ we get the following useful equation with 
details presented in the appendix 


H 0 D l 


c(l + z) . [ y/\^k\ f Zl 

•\Z|Hfc| y/^Kn Jo 


— ==dz (16) 

+ ^ + + £d + + 


and sinn is sin for < 0 and sinh for > 0. Note that within the 
denominator of the integral Ha appears in only linear combination with z 
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as opposed to appearance of this term in the more conventional use (Carroll 
et al. 1992, Tonry et al., 2003) in this equation 


H 0 D l 


c( 1 + z) 
\/|^fc| 


sinn{|rifc| 


1/2 



1 

\/(l + z) 2 (l + £l m z) 


= =dz. 

z(2 T z)H A 


(17) 

Again, the nature of sinn depends upon spacetime curvature as above and 
equations (16) and (17) are different forms of the same solution. 


4 RESULTS 

There are some problems of inconsistency that arise with this solution of 
the standard model and we observe these near parameter values presently 
considered important. We find the root portion of the denominator of (10) 
causes these inconsistencies across a range of normalized matter densities 
as displayed in Figure 1 for Q m from 0.30 to 0.05. These graphs present 
the discontinuity, with respect to equation (10), increasing in intensity with 
corresponding value of z as il rn approaches 0. For instance, at f l m of 0.30, 
currently of interest (Feldman et al., 2003), fi A appears inconsistent from 
1.713 to 1.714 with the associated values of H*. of -1.013 and -1.014. This 
inconsistency widens considerably at H m of 0.05, perhaps of interest (Fall, 
1975), where at z ~ 2.65 the inconsistency in fl A lies from 1.224 to 1.250; 
though the breadth of the inconsistency increases with increasing z and 
decreasing H A the exact values of these inconsistencies decrease along with 
disappearing matter density. 

We have traced these inconsistencies at H m of 0.05 and 0.30, with various 
values for H A and present results in Figures 2 and 3 as plots of HqDl with 
respect to log z. Both figures demonstrate for large values of the cosmologi¬ 
cal constant, 1.2 to 1.7, and with reasonable values of spacetime curvature 
Hfc drastic deviations from useful solutions of HqDl. At the lower matter 
density of 0.05 (Fig. 2) this function does not straightforwardly correlate 
HqDl with z, for values of H A above 1.2 and for z above about 2. Here the 
critical value for H A is a rather low 1.224. With increasing H A from 0.70 to 
0.95, we find the larger values of z correlate with the smaller values of H A 
as should be observed in a Universes with increasing matter densities and 
decreasing anti-matter densities. 

As plotted in Fig. 3 (H m = 0.30), all curves are shifted to higher values of 
z, as expected for a Universe more populated with matter than in Figure 
2. The model appears well behaved with larger values of H d Dl correlating 
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Redshift z Redshift z 


Figure 1: {U m (l + z) 3 + fl A + f4(l + z) 2 } versus redshift z for various 
normalized matter densities between 0.30 and 0.05. 

with larger values of redshift up to perhaps Q A of 1.2; solutions of H q Dl at 
values of fi A greater than this and above z of 3.25 are questionable. Here, 
evaluations of HqDl with smaller values of Q A appear indiscriminate but 
perhaps useful up to z of about 3, but no further. Beyond this redshift z is 
predicted to correlate with anti-gravity density. Though these problems of 
inconsistency are for f2 A greater than 0.70 - just outside the range currently 
thought to be important - as the Universe expands, f2 A is supposed to dom¬ 
inate with f i m eventually approaching 0 and the standard model may fail 
for even moderate values of the redshift and Q A . 

We have found the FRW model does not obviously fail for negative values 
of U A over a much wider range of Q m and z. Figure 4 presents log HqD ^ 
versus logz for Q m of 0.05 with U A from +1 to -2. As expected for a Universe 
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Figure 2: Plots of HqDl as functions of log z with various values of H A with 
Q m of 0.05 and c normalized to 1. 

dominated by attractive forces, HqDl decreases with decreasing f1 A and z 
decreases with decreasing Q A in smooth manners because the Hubble flow 
must increase with increasing drag on expansion. Note that a positive cos¬ 
mological constant (fi A = 1) does follow this regular pattern above redshift 
of 6.59 in the Universe of low matter density (0.05); greater than this dis¬ 
tance the energy is predicted to become strongly attractive. This ordering 
also exists at of 0.30 (Fig. 5) but is more pronounced; up to redshifts of 
about 3.25 a U A of +1.0 acts as anti-gravity, but at greater distance than 
this fl A acts attractively. 

5 DISCUSSION 

The selection of a solution to the standard model of cosmology using an el¬ 
liptical form allows us to probe the ability of the model to estimate Hubble 





Figure 3: Plots of HqDl as functions of log z with various values of II A with 
n m of 0.30 and c normalized to 1. 


expansion at great distances with dependence upon fi A . We have found, for 
values of J) A between 0 and -2 and values of Q m currently thought realistic, 
the standard model returns a consistent value for HqDl at z < 20. This 
is greater than redshifts currently reported for the farthest galaxies, though 
this value may be surpassed soon. When this model is solved for values of 
0 A greater than 0 problems arise at high values of 2 , where positive J) A acts 
as attractive. In regions of fI A greater than +1 inconsistencies are observed 
of a general nature which limits the usefulness of the standard model to 
regions of high matter density - densities consistent with the Einstein-de 
Sitter model and currently thought important. If the Universe is rushing 
towards a state of low matter density, increasing dependence upon fi A with 
a trace of closed curvature (preferred by Einstein among others) the stan¬ 
dard model will probably fail. These regions of inconsistencies might be 
interpreted by some, as evidence for an epoch of matter dominated Universe 
with a low Universe expansion rate followed by the current epoch of vacuum 
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Figure 4: Log HqDl as a function of log z for various values of 12 A from 
+1.0 to -2.0 with 12 m of 0.05 and c normalized to 1. Note the ’’cross-over” 
of the trace of the solution with Q A of +1.0. 


energy domination and faster expansion, though it seems to us that earlier 
epochs might be best not judged using the standard model. Very unfortu¬ 
nately, this also means calculations of epochs of newly discovered galaxies 
exhibiting high redshifts will remain unreliable until the FRW model can be 
better adapted for greater distances. The standard model has limitations 
of usefulness with respect to accurate predictions of large galactic distances 
and Universe age. Inclusion of positive values for the vacuum energy further 
restricts the range of useful predictions. Other models or model variations 
incorporating modified application of the vacuum energy, but not suffering 
regions of inconsistency, might be preferred. We hope to introduce one such 
model soon. 
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Figure 5: Log HqDl as a function of log z for various values of 17 a from 
+1.0 to -2.0 with f7 m of 0.30 and c normalized to 1. Note the ’’cross-over” 
of the trace of the solution with of +1.0. 
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APPENDIX 


We derive equation (16), the elliptical form useful for astronomy from our 
equation (9). The FRW metric, allowing a(t) = -g- is 


dR 0 r 

dt 


R °dt ~ 


— (1 — Ax 2 ) 1 / 2 

R 

-L(l - kr^ 
a (t) 


with rearrangement becomes 


R 0 dr dt 

(1 — kr 2 ) 1 / 2 a(t) 


(18) 

(19) 


( 20 ) 


Remembering that r= H a t we multiply through by H 0 to get 


H 0 R 0 


dr 

(1 — kr 2 ) 1 / 2 


dr 

a(t) 


( 21 ) 


and multiply through again using the relationship for the normalized Pj. = 

k 


'WR 


H n R n 


dr 


dr 


\l + n o R 2 H 2 r 2 y/ 2 a(t) 
Substituting equation (9) for dr we have 


( 22 ) 


H 0 R 0 


dr 

(1 + n o R 2 H 2 r 2 y/ 2 


1 y/ a{r)da{r) 

a(r) y/Vt m + P A a(r) 3 + Pfc«('r) 


(23) 


and integrating both sides and rearranging some constant parameters 


H 0 R 0 


r r i 


1 


. /i i o r?2 


-dr = 


(24) 

then changing the variable on the left side using y = y/TikR 0 H 0 r and using 
the redshift relationship a = ■prprj- for substitution of the right hand side we 
get 


r^/n^Ro Hqti 


\Z^k ■ 


V 1 + y 2 


dy = - 




Z1 (! + Z) 2 y/P m + 


dz. 


(25) 
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The integral on the left hand side is arcsinh(y) so the above becomes 


— arcsin h(^yn^R 0 H 0 n) 

v“fc 



1 

y / fi m (l + z) 3 + + f4(l + z) 2 


dz 


(26) 

and substituting to arrive at the measurable Dl using the relationships 
D m = R 0 r\ and Dl = (1 + z)D m we arrive at 


H 0 D l 


i+z . fV w r 

— ; - sin nl . / 

\/|14| VTU J o 


V /( l +-*) 3 + 


1 




+ £*-(1 +zf 


-.dz 


(27) 


and with introduction of c for the speed of light becomes our equation (16). 


14 



